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Abstract 

O 

C 

We present multi-model global datasets of nitrogen and sulfate deposition covering 
time periods from 1850 to 2100, calculated within the Atmospheric Chemistry and 
Climate Model Intercomparison Project (ACCMIP). The computed deposition fluxes 
5 are compared to surface wet deposition and ice-core measurements. We use a new 
dataset of wet deposition for 2000-2002 based on critical assessment of the qual- 
ity of existing regional network data. We show that for present-day (year 2000 AC- 
CMIP time-slice), the ACCMIP results perform similarly to previously published multi- 
model assessments. For this time slice, we find a multi-model mean deposition of 
io 50Tg(N)yr _1 from nitrogen oxide emissions, 60Tg(N)yr _1 from ammonia emissions, 
and 83Tg(S)yr _1 from sulfur emissions. The analysis of changes between 1980 and 
2000 indicates significant differences between model and measurements over the 
United States but less so over Europe. This difference points towards misrepresenta- 
tion of 1980 NH 3 emissions over North America. Based on ice-core records, the 1850 
15 deposition fluxes agree well with Greenland ice cores but the change between 1850 
and 2000 seems to be overestimated in the Northern Hemisphere for both nitrogen 
and sulfur species. Using the Representative Concentration Pathways to define the 
projected climate and atmospheric chemistry related emissions and concentrations, 
we find large regional nitrogen deposition increases in 2100 in Latin America, Africa 
20 and parts of Asia under some of the scenarios considered. Increases in South Asia 
are especially large, and are seen in all scenarios, with 2100 values more than double 
2000 in some scenarios and reaching > 1300mg(N)m _2 yr 1 averaged over regional 
to continental scale regions in RCP 2.6 and 8.5, ~30-50% larger than the values in 
any region currently (2000). The new ACCMIP deposition dataset provides novel, con- 
25 sistent and evaluated global gridded deposition fields for use in a wide range of climate 
and ecological studies. 
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1 Introduction 


g 
(/) 
O 
c 

The global nitrogen cycle is of importance for a number of key-issues, such as ecology 
and biodiversity (e.g. Phoenix et al., 2006; Bobbink et al., 2010; Butchart et al., 2010), 
eutrophication and acidification (e.g. Bouwman et al., 2002; Rodhe et al., 2002; Fisher 
5 et al., 2011), climate change-carbon cycle interactions (e.g. Thornton et al., 2007; Reay 
et al., 2008; Zaehle et al. 2010), food and energy production (Sutton et al., 2011). 
Nitrogen emissions also impact human health through particulate matter and ozone 
formation. Clearly, nitrogen is central to many aspects of life on Earth (Galloway et al., 

2008; Fowler et al., 2012). 

io Dinitrogen (N 2 ) is the most abundant component in the atmosphere, but is chemically 
unreactive, only contributing to tropospheric chemistry through lightning and associ- 
ated NO formation (Franzblau and Popp, 1989). More reactive atmospheric nitrogen 
components include oxidized NO y (= NO + N0 2 + other minor inorganic components 
and organic nitrogens), the long-lived greenhouse gas dinitrous oxide (N 2 0; 131 ±10yr, 

15 Prather et al., 2012), and reduced nitrogen NH X (= NH 3 + NH 4 ) as well as organic ni- 
trogen components such as amines (Kanadidou et al., 2012). NO y and NH X are collec- 
tively identified as reactive nitrogen (N r ). Anthropogenic emissions of N r components 
are estimated (van Aardenne et al., 2001 ; Lamarque et al., 201 0) to have increased by 
a factor of 5 to 10 since 1850. Knowledge on the link between N r generation by human 
20 activities and its subsequent impacts requires an accurate description of largescale 
emissions, atmospheric chemistry transport and removal, all of which occur on rela- 
tively fast timescales compared to nitrogen in the other compartments of relevance to 
the Earth system: terrestrial, coastal zones and open ocean. 

Global atmospheric-chemistry transport models are routinely used as a tool to calcu- 
25 late the global dispersion of N r . Since individual models are prone to specific errors and 
multi-model means typically outperform individual models (e.g. Reichler and Kim, 2008; 
Shindell et al., 2012), it has become common to use model ensemble calculations to 
improve the quality of the calculations. Previous deposition ensemble studies include 
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Lamarque et al. (2005) using 6 models and focusing on NO y deposition, Dentener et al. 
(2006) using 23 models and discussing NO y , NH X and SO x deposition under current 
and 2 future scenarios, and Sanderson et al. (2008) using 15 models and focusing on 
export of NO y and subsequent deposition. 

5 The present study uses deposition fields generated by 10 models that participated 
in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP; 
Lamarque et al., 2012). ACCMIP was designed to inform aspects of the forthcom- 
ing Intergovernmental Panel on Climate Change Assessment Report #5 regarding 
the role of long-term changes in atmospheric chemistry on climate and vice-versa. 
io In this ACCMIP study simulations were provided for time slices representative of the 
period around 1850, 1980, 2000, 2050 and 2100 (see Sect. 2), where the future atmo- 
spheric chemistry-climate conditions follow the 4 Representative Concentration Path- 
ways (Moss et al., 201 0). Because of the importance of sulfur deposition (in association 
with nitrogen deposition) for the understanding of soil and water acidification (Fisher 
15 et al., 2011), we have included an analysis of sulfur deposition to the more traditional 
nitrogen deposition analysis. The models used in this analysis represent a combina- 
tion of the current generation of Chemistry-Transport Models and Chemistry-Climate 
Models, with somewhat refined horizontal and vertical resolution and more detailed 
descriptions of chemical processes compared to previous studies (namely the results 
20 discussed in Dentener et al., 2006; Sanderson et al., 2008). Rather than represent- 
ing specific meteorological years - as in previous studies - each model meteorology 
is rather representative of the decade under consideration (e.g. 1850-1859; or 2000- 
2009) and hence include the effect of climate change. 

This paper is structured as follows. Section 2 gives a short description of the models 
25 and the ACCMIP experiments. We also define there the methodology for computing 
the multi-model mean (MMM). Section 3 describes the measurement datasets used 
in this paper, from surface sites and from ice cores. In Sect. 4, we evaluate the multi- 
model mean for the 2000 period - when extensive measured deposition datasets are 
available, and which allows us also to make a retrospective analysis of the quality 
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of data in this model dataset compared to the earlier and widely used Photocomp 
(Dentener et al., 2006) and HTAP (Sanderson et al., 2008; Vet et al., 2013) deposi- 
tion datasets. In Sect. 5, a more limited analysis is given for the changes in deposi- 
tion (and its drivers) since the 1980s and, using ice-cores, for the change between 
decades centered around 1850 and 2000. The 1980-2000 period is of particular inter- 
est since, due to implementation of national and international emission control mea- 
sures, strong emission reductions have been reported, especially over North America 
and Europe, with commensurate consequences for deposition. Section 6 describes the 
overall global structure of total deposition (dry and wet) from 1850 to 2100. A discus- 
sion and conclusions follow in Sect. 7. 

c n 

O' 

13 

2 ACCMIP simulations and multi-model mean 

■o 

CD 

—i 

An overview of the Atmospheric Chemistry and Climate Model Intercomparison Project 
(ACCMIP) simulations and participating models is given in Lamarque et al. (2012). 
Consequently, we focus here on the aspects relevant to the analysis of nitrogen and 
sulfate deposition. 

ACCMIP provides analysis on the role of atmospheric chemistry changes in the near- 
term-climate forcing included in the CMIP5 (Climate Model Intercomparison Project 
Phase 5; Taylor et al., 2012) simulations, including the chemical composition changes 
associated with the CMIP5 prescribed forcings. The ACCMIP simulations used in the 
present study (Table 1) consist of time slice experiments (for specific periods span- 
ning 1850 to 2100 with a minimum increment of 1 0 yr) with chemistry diagnostics. 

Each requested simulation is labeled as Primary (“P”) or Optional (“O”). Although sim- 
ulations using the Representative Concentration Pathway 6.0 (RCP6.0, Masui et al., 

201 1 ) emissions were performed, not enough groups provided the necessary fields for 
a meaningful analysis of nitrogen and sulfur depositions; therefore, this work focuses 
on the remaining three RCP projections RCP2.6, RCP4.5 and RCP8.5 (see van Vu- 
uren et al., 2011 and references therein). A primary output of these simulations was 
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nitrogen and sulfate (dry and wet) deposition fields (see Table SI in Lamarque et al., 
2012). Although the ACCMIP models were required to specify the anthropogenic and 
biomass burning emissions according to Lamarque et al. (2010) for 1850-2000 and 
the RCP emissions (van Vuuren et al., 2011) beyond, there is a range of emissions 
that were used in ACCMIP models, mostly from variations in the treatment of natural 
emissions (Lamarque et al., 2012; Young et al., 2013). 

As a first step in the quality control of the model calculated depositions, we ana- 
lyze for each model the balance between emission and deposition as reported by the 
modeling groups (Tables 2-4) for the year 2000 time slice experiment. This analysis 
is done separately for nitrogen originating from nitrogen oxide emissions (deposition 
fields are referred as a group as NO y ) or ammonia emissions (NH X ) and sulfur emis- 
sions (SO x ). We find that the NO v deposition is larger than the emissions (surface and 

y -i “ 

upper-air, including lightning) by approximately 1 Tg(N)yr , representing the input of 
nitrogen (mostly in the form of nitric acid) from the stratosphere, in general agreement 
with observational estimates of 0.45Tg(N)yr _1 from Murphy and Fahey (1998), except 
for the GISS model where this influx is approximately 5Tg(N)yr _1 . In the case of am- 
monium, the balance between emissions and deposition is clearly attained. In the case 
of sulfur, this balance cannot be fully evaluated due to the lack of diagnostics on the 
formation of sulfate aerosols from the emitted dimethylsulfide (DMS). Boucher et al. 
(2003) estimated the yield (DMS to sulfate aerosols) to be 87% when both gas-phase 
and aqueous-phase reactions are taken into account, but this number will be some- 
what model dependent. Within that limitation, it is reasonable to assume that balance 
between sulfur deposition and emission is achieved for the listed models. 

The focus of this paper is on documenting the multi-model mean (MMM) generated 
for each time slice. The MMM is constructed by linearly interpolating the model gen- 
erated monthly fields (for example wet deposition combined for all NO y species) at 
their native horizontal resolution (Table SI ) to a common 0.5° x 0.5° grid (finer than any 
model grid), identical to the emission grid used in Lamarque et al. (2010). Then, each 
field is averaged across models at the original monthly temporal resolution to generate 
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its multi-model mean and standard deviation (Tables S2-S4). As indicated in these ta- 
bles, the largest number of models (up to 9) generated the necessary fields is found 
for the historical NO y deposition, followed by the RCP8.5 simulations. The number of 
models performing ammonium chemistry and deposition is however much smaller (2- 
3 models, depending on the simulated time period) while between 2 and 6 models 
have provided sulfate fields. In all cases, the MMM is constructed using all available 
model results, therefore leading to variations between the specific models used in the 
average. 

</)’ 

O 

c 

3 Description of observational deposition data 

13 

Under the auspices of the World Meteorological Organization (WMO) a precipitation 
chemistry expert group has performed a critical analysis of available wet deposition 
data for the years 2000 to 2002. While dry deposition may be available at specific sites, 
this is not a directly observed quantity and we therefore do not use this information. We 
use the WMO-processed wet deposition datasets in our analysis of the performance of 
the MMM in the 2000 time slice. 

In the WMO dataset, a careful analysis is made of worldwide reported data of 
wet precipitation chemistry, and data quality qualifiers are provided to deposition 
data obtained mainly from networks in Europe European Monitoring and Evaluation 
Programme, EMEP; http://www.nilu.no/projects/ccc/emepdata.html; the United States 
(National Atmospheric Deposition Program, NADP; http://nadp.sws.uiuc.edu/NTN/ 
ntnData.aspx); Canada (Canadian Air and Precipitation Monitoring Network, CAP- 
MoN; http://www.ec.gc.ca/rs-mn/default.asp?lang=En&n=752CE271-1); Asia (mainly 
fromtheAcid Deposition Monitoring Network in East Asia EANET; http://www.eanet.cc), 
and Africa (Deposition of Biogeochemically Important Trace Species (DEBITS), http: 
//debits. sedoo.fr). We will focus our analysis on the NADP, EMEP and EANET obser- 
vations. For more detailed information on data networks, data selection criteria, and 
an application to the HTAP Phase 1 deposition dataset, we refer to Vet et al. (2013), 
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as well as the data networks indicated above. The evaluated wet deposition datasets, 
soon to be available through the World Data Centre for Precipitation Chemistry, only 
contain a subset of the available stations for that period, i.e. stations that correspond to 
regionally representative sites that fulfilled specific data completeness and quality con- 
5 trol measures. In this paper, we use the wet deposition measurements of nitrate (NOg), 
ammonium (NH 4 ) and non-sea salt sulfate (SO^), which were qualified as “good”. 

Over Europe (EMEP) and the United States (NADP), deposition measurement sites 
typically started around the 1980s. We use these early measurements to evaluate the 
change in deposition over these two decades. To this purpose, we computed from the 
10 raw data (available at the aforementioned web sites), and using the filtering protocols 
defined by the specific networks, two sets of 6 -yr averages (1980-1985 and 1997- 
2002 ) for sites that had sufficient observations for both time slices. 

Although a quality control similar to the WMO evaluation was not available for the 
1980s data, comparison of the WMO processed 2000-2002 data with our 1997-2002 
15 indicate very good agreement over the respective regions (not shown), validating the 
suitability of our processed dataset to study trends in deposition from the 1980s to the 
2000 s. 

C 

Historical records of nitrate, ammonium, and sulfur were also developed from high- 
depth-resolution measurements in ice cores (see Table S5 for their geographical infor- 
20 mation) using an established continuous ice core analytical system (e.g., McConnell 
and Edwards, 2008). Nitrate and ammonium were measured using spectrophotome- 
try and fluorimetry, respectively, with standard flow through methods (Roethlisberger 
et al., 2000). Total sulfur concentrations were measured using magnetic sector Induc- 
tively Coupled Plasma Mass Spectrometry (McConnell and Edwards, 2008). At core 
25 sites with sufficient annual snowfall, the ice core records were dated using annual layer 
counting of a range of seasonally varying chemical species (Sigl et al., 2012). Synchro- 
nization to well-known volcanic layers was used for dating at core sites with very low 
snowfall (e.g., East Antarctica) where annual layers are not preserved (Anschutz et al. 
2011) or at lower elevation sites (e.g., Akademmi Nauk, McCall Glacier, Flade Isblink) 
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where surface melting and percolation make annual layer identification difficult (Opel 
et al., 2009). Uncertainty (la) in the dating is estimated at ±1 yr for the annually dated 
sites and ±3yr for all other sites. Decadal averages centered on 1855 and 2000 were 
computed from the high-resolution measurements. 

"O 

CD 

— t 

4 Evaluation of ACCMIP year 2000 time slice deposition 

While observation and model data are available at the monthly time resolution, we fo- 
cus our analysis on the annual mean. This choice is made to limit the discussion to the 
long-term trends in nitrogen and sulfate deposition. In addition to the present intercom- 
parison project, MMM results from two previous studies are available for comparison 
and analysis of potential improvements: PhotoComp (Dentener et al., 2006) and HTAP 
(Sanderson et al., 2012; Vet et al., 2013). These previous studies are partially inde- 
pendent of the present one, with different emissions and a different sets of models or 
different versions of the same models. Using the WMO dataset and the model results 
interpolated to the location of the observing stations, we can statistically analyze the 
ability of the models to reproduce the observations. 

We display in Fig. 1 the global and regional distributions of N0 3 NOg (Fig. la), 

NH 4 (Fig. 1b) and SO^ - (Fig. 1c) wet deposition in the MMM compared to the WMO 
dataset (over Europe and the United States). Deposition is clearly strongly correlated 
with emissions (i.e. largest in the Northern Hemisphere, but also with larger amounts in 
areas of biomass burning and large soil emissions such as Central Africa), albeit with 
significant downwind propagation for each compound. Unlike the nitrogen sources, 
sulfur emissions from degassing volcanoes (Andres and Kasgnoc, 1998) lead to the 
formation of deposition hotspots in areas with low anthropogenic emission levels such 
as Central America. 

In the case of NOg deposition, we find that all 3 experiments (ACCMIP, HTAP and 
PhotoComp) perform rather similarly (Fig. 2 and Table 5). However, none of the model 
simulations are able to capture the observed high deposition rates over East Asia or 
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Europe, possibly due to the proximity of the observing stations to significant sources, 
features that cannot easily be captured with the coarse grid of presently used models 
(Zhang et al., 2012). In contrast, in the high emission region of Northeastern United 
States, there are regions of overestimation as well as underestimation. 

In the case of NH 4 , the ACCMIP results are slightly lower than the previous stud- 
ies over both the North American NADP and European EMEP domains and there- 
fore with a larger negative bias. Because of errors in the regridding of the NH 3 agri- 
cultural emissions over China in ACCMIP (see http://www.iiasa.ac.at/web-apps/tnt/ 
RcpDb/dsd?Action=htmlpage&page=download), the performance in ACCMIP (Fig. 2) 
is considerably worse than the previous studies over Asia, with a significant underes- 
timate (by approximately 30%, see Table 5) for the highest deposition rates over East 
Asia. 

The sulfate (SO^ - ) deposition is somewhat improved in ACCMIP over the NADP sta- 
tions, while the very high biases previously found in the HTAP dataset over the EMEP 
domain are now considerably reduced, so that the ACCMIP models might actually be 
underestimating wet deposition over that region. The deposition over EANET is char- 
acterized by a larger negative bias in ACCMIP than PhotoComp or HTAP, although the 
overall correlation remains high. 

Overall, the performance of the ACCMIP MMM is similar to PhotoComp and HTAP, 
suggesting that no significant improvement or worsening has been made since those 
2 studies in the representation of emissions, chemical processing and deposition pro- 
cesses. The analysis also shows that ammonium deposition over East Asia is most 
likely underestimated. On the other hand, for all considered species, the ACCMIP MMM 
deposition tends to be slightly overestimated over the NADP stations. 

U) 

O 

c 
( / > 

5 Evaluation of historical deposition changes: 1980-2000 and 1850-2000 

Tl 

To identify potential changes in deposition over time, we use the 1980 and 2000 time 
slice experiments. Unfortunately, neither HTAP nor PhotoComp provide historical time 
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slices and therefore such changes cannot be discussed within that context. Over that 
period, Europe and North America have seen significant changes in nitrogen and sulfur 
emissions, from a combination of changes in anthropogenic activities and air quality 
regulations (Granier et al., 2011; Xing et al., 2012). 

5 To limit issues associated with interannual variability and uneven time sampling from 
the observed wet deposition rates, each station much have at least 36 months of avail- 
able data for each 6-yr averaging period (see Sect. 3). 

Over the United States, the NADP measured nitrate deposition change between 
1980 and 2000 consists mostly of reductions, especially in the eastern portion of the 
io United States (Fig. 3a). However, the MMM shows only a minor change over the west- 
ern United States, and actually simulates a minor increase over the eastern States. In- 
deed, all models except GISS (Table 6) indicate a relative increase in nitrate deposition 
averaged over the United States. This simulated increase in wet deposition is a com- 
bination of minor increases in NO x emissions (Fig. SI a) and in precipitation amounts 
is (Table 6). Based on the CPC Merged Analysis of Precipitation (CMAP; Xie and Arkin, 
1997), there is however indication of a small observed decrease in precipitation be- 
tween 1980 and 2000 over the United States (Fig. S2), while the MMM has limited 
inter-model agreement and actually shows a slight increase in precipitation when av- 
eraged over the whole United States (Table 6). It must however be recognized that 
20 the changes over the United States are relatively small and may therefore be strongly 
affected by interannual variability. Over the EMEP network, many sites indicate a very 
strong reduction consistent with emission change (Fig. SI a), with the exception of for- 
mer Yugoslavia and the southern tip of Norway. The MMM captures well the strong 
reduction, with however a smaller amplitude than is observed at most locations. In this 
25 case (Table 6), all models agree on the change in emissions (-12% for the mean) and 
all but MOCAGE show an increase in precipitation. 

Measured ammonium deposition (Fig. 3b) changes over the NADP network show 
a mixture of increases and decreases, with again the largest decreases over the North- 
eastern United States. Similar to nitrate deposition, the MMM is much more uniform 
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than observed and only shows a small decrease in the mid-Western United States. 
This is due to the fact that the anthropogenic emissions of NH 3 are almost identical be- 
tween 1980 and 2000 (see Fig. Sib). This is possibly erroneous but more analysis is 
beyond the scope of this paper. Over Europe, the EMEP station data indicate an overall 
5 strong decrease, with the exception of sites in France, Italy and Norway. In contrast, 
the MMM indicates an overall small increase over Western Europe and a decrease 
over Eastern Europe, consistent with the significant emission change (Fig. Sib). Local 
factors, not included in the models, play a role in these differences. 

Over the United States, except for a dozen sites scattered over the NADP network, 
io observed sulfate deposition east of 100°W (Fig. 3c) is characterized by large reduc- 
tions driven by emission change (Fig. Sic). The MMM is however underestimating the 
amplitude of the changes, especially for the largest changes. Over the EMEP network, 
the MMM captures well the general deposition decrease including its largest change 
over Germany, clearly in agreement with the emission change (Fig. Sic), 
is Overall, the analysis of changes between 1980 and 2000 shows limited agreement 
between MMM and observations over the United States. In particular, the observed 
changes over the NADP network indicate much higher variability and amplitude than 
is simulated, driven by fairly small changes in emissions. Also, this analysis highlights 
the potential error in the NH 3 emission change over the United States between 1980 
20 and 2000. While the deposition change over Europe from the MMM is not as large as 
observed, the overall patterns of change are much better represented, except for NH X 
deposition over Western Europe. This suggests that NO x (and to a lesser extent SO x 
and NH X ) emissions changes may better captured over Europe in the historical emis- 
sion dataset used in ACCMIP (Lamarque et al., 2010), since emission change is the 
25 main driver for deposition change. However, further analysis is required for a complete 
understanding of the applicability of deposition data to constrain emission inventories. 

In addition, we use ice core records of deposition over the Northern Hemisphere and 
over Antarctica to study the ability of the models to represent changes since 1850. To 
limit the importance of interannual variability in the ice-core record, we use the average 

6259 


ACPD 

13, 6247-6294, 2013 


Multi-model mean 
nitrogen and sulfur 
deposition 

J.-F. Lamarque et al. 


Title Page 

Abstract I Introduction 
Conclusions I References 
Tables I Figures 



Printer-friendly Version 
Interactive Discussion 




















values for 1850-1860 and 1995-2005 to compare against the 1850 and 2000 time 
slices, respectively. As a basic test of the applicability of the MMM to the polar regions, 
we show in Fig. S3 a comparison of year 2000 precipitation against the gauge-based 
climatology of Yang et al. (2005). We find that the IVIMM compares rather well to the 
5 observations throughout the Arctic circle, with the caveat that no rain gauge is available 
over Central Greenland. Additional analysis of precipitation over the Arctic is provided 
in Lee et al. (2013). No equivalent data are available for Antarctica. 

We find that there is a strong agreement between the 1850 observed and simulated 
nitrate depositions (Fig. 4a) for 1850, both in regional structure and intensity. Nitrate 
io deposition in 2000 is however overestimated at all Northern Flemisphere sites but one. 
On the other hand, the year 2000 simulated deposition over Antarctica, which is much 
less affected by changes in anthropogenic NO x emissions but still shows an overall 
increase in nitrate deposition, is in good agreement with the observations. 

Ammonium deposition (Fig. 4b) is strongly overestimated at the Greenland sites in 
15 both 1850 and 2000 time slices (and so is the change 2000 minus 1850; not shown). 
Wolff et al. (201 3) argue that it is actually difficult to derive a change of measured NFI 4 
deposition over Greenland, due to the large dependency on the highly variable contri- 
bution of biomass burning to NFI 4 deposition. On the other hand the other two sites 
in the Northern Flemisphere (McCall, Alaska and Akademmi Nauk, Siberia; see Table 
20 S5 for their exact locations) are well captured. This is true for the 2000 time slice as 
well, where the McCall glacier shows a decrease since 1850, possibly associated with 
a change in biomass burning emissions over North America. Over Antarctica, ammo- 
nium deposition rates over the Western portion of the ice shelf are well represented, 
while they are somewhat overestimated in the eastern sector. This is the case for both 
25 1850 and 2000. Increases in the simulated deposition between 1850 and 2000 seem 
to be larger than observed, albeit the levels of deposition are small compared to those 
observed in the Northern Flemisphere. 

Sulfate deposition in 1850 (Fig. 4c) is represented quite well at the Northern Flemi- 
sphere and Antarctic sites, except for McCall, Alaska and central Greenland (Tunu 
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glacier). Sulfate deposition in 2000 is overestimated at all ice-core sites of the North- 
ern Hemisphere. Sulfate deposition over Antarctica has a well-defined east-west sep- 
aration, well captured by the model. There is little variation between 1850 and 2000 
since the sulfate production is primarily driven by DMS oxidation. In many models used 
in ACCMIP, the DMS source does not vary with time, although Cameron-Smith et al. 
(2011) suggest it should, and changes in DMS emissions could represent an important 
feedback in the climate system (Charlson et al., 1987). 

Overall, the ice-core comparison of the MMM indicates a reasonable representation 
of the pre-industrial (1 850) conditions but tends to overestimate the present-day (2000) 
conditions. Similar conclusions are found for 1980 (not shown), indicating that biases 
in the year 2000 deposition fields are not related to transient features in ice formation. 

"O 

0 ) 

T3 

6 Total nitrogen and sulfate deposition 

In this section, we document the MMM global and regional distributions of total deposi- 
tion (wet + dry) of NO y , NH X and SO x and their changes from 1 850 to 21 00 under the 
RCP2.6, 4.5 and 8.5 projection scenarios (Fig. 5 and Tables S2-S4). 

The historical increase of NO y deposition (Fig. 5a) since 1850 took place mostly 
in the Northern Hemisphere. It is characterized by deposition rates larger than 
1 kg(N) ha -1 yr -1 over most of the continental areas, the North Atlantic and the out- 
flow oceanic areas of Eastern Asia, the Indian subcontinent and Central Africa. By 
2030, the projections over the United States and Europe are quite similar for RCP2.6 
and RCP4.5, while RCP8.5 shows a smaller reduction from the 2000 levels. In East- 
ern Asia, deposition of NO y is showing levels above 2000, with the largest increase 
seen in RCP8.5. On the other hand, over the Indian sub-continent the largest depo- 
sition increases were found in RCP4.5. Only small changes occur in the Southern 
Hemisphere. By 2100, most of the continental areas, except for India, are projected to 
return to pre-1980 levels of deposition. Another exception is the larger deposition over 
the Northern Pacific ocean in RCP8.5, consistent with NO x emissions from shipping. 
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Since the MMM NO y deposition is computed from the largest number of models, it is 
reasonable to document intermodel variability for that diagnostic (Fig. S6). We find that 
over most continental areas and for most time slices, the inter-model standard devi- 
ation is in the 10-20% range, with Central Asia and South America being somewhat 
5 larger (20-30%). 

While the increases in nitrogen deposition from ammonia emissions between 1850 
and 2000 affect the various regions similarly to the impact of NO x emissions (Fig. 5b), 
the RCP projections for NH 3 emissions are indicative of a very different trajectory 
(Lamarque et al., 2011), with continuous increases over most continental and oceanic 
io regions, with the exception of the oceans south of 30° S. These deposition changes are 
mainly driven by the projected increases in inorganic fertilizer use needed to provide 
more and higher quality food for a growing worldwide population, with no policies in 
places to abate the emissions. 

Simulated sulfate deposition reached higher levels and was more widespread in Eu- 
15 rope than North America in 1980, but was similar in both regions by 2000. The only 
significant increase between 1980 and 2000 is over China, associated with its increas- 
ing use of coal for power generation (Smith et al., 2011). The projected changes over 
these regions indicate a gradual phase-out of anthropogenic emissions by 2100, in- 
cluding over China. Even by 2030, emissions over China are projected to be no larger 
20 than in 2000, a trend that might be reflected in the most recent estimates over China 
(Smith et al., 2011; Klimont et al., 2013). Similar to NO y deposition, only the Indian 
subcontinent is projected in RCP8.5 to have sulfate deposition levels higher in 2100 
than in 1980. Sulfate deposition over the oceans is also considerably reduced in the 
RCP scenarios, but depends on the specific scenarios for their projection of shipping 
25 emissions. 

The combination of deposition fluxes as 2-SO x + NO y -NH x provides an indication 
of the degree of acidity (Fisher et al., 2011) contained in the deposited fluxes (Fig. 6). 
In particular, it clearly shows that, over the 21st century, continental areas over the 
Northern Hemisphere will have a tendency towards more basic deposition with the 
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increase NH 3 emissions. This analysis also clearly shows the spurious NH 3 emissions 
over Mongolia (as discussed in Sect. 4). 

Following the discussion of averaged deposition rates over various regions in Lamar- 
que et al. (2011 see their Tables 5 and 6; note that the NCAR-CAM3.5 results used in 
5 the present ACCMIP analysis are equivalent to results in those), we present the same 
diagnostics here (in the same units), with the exception of RCP6.0 (Tables 7 and 8). 
Overall, the total (from NO y and NH X combined, since both affect vegetation in the 
same fashion) nitrogen deposition is expected to increase between 2000 and 2100 
over many regions, especially in the case of RCP8.5. Furthermore, only Europe and 
io North America are projected to see a reduction in their nitrogen deposition. Sulfate 
deposition rates are clearly projected to significantly decrease by 2100, as expected 
from Fig. 5. All these confirm the results published in Lamarque et al. (201 1 ) in a single 
model study using the same emission fields applied in the present study. 

7 Discussion and conclusions 

g 

</)' 

is We have presented in this paper the multi-model mean annual deposition of nitrogen 
and sulfate as simulated in the Atmospheric Chemistry and Climate Model Intercom- 
parison Project (Lamarque et al., 2012). We have made considerable use of network- 
based wet deposition datasets (including a new dataset based on expert analysis of 
the deposition data; Vet et al., 2013) and ice-core records to evaluate deposition rates 
20 and changes in the deposition since 1980 and 1850. We have found that the ACCMIP 
multi-model mean behaves rather similarly to previous multi-model analysis (Dentener 
et al., 2006; Sanderson et al., 2008), with the notable exception of ammonium depo- 
sition over Asia, which is considerably worse in ACCMIP due a regridding error in the 
underlying emissions over that region (Lamarque et al., 2010). 

25 Beyond the present-day analysis, we discuss a comparison of the change in depo- 
sition rates between 1980 and 2000 (using surface wet deposition) and between 1850 
and 2000 (using ice-cores). Although the deposition in 2000 is rather well simulated in 

6263 


ACPD 

13, 6247-6294, 2013 


Multi-model mean 
nitrogen and sulfur 
deposition 

J.-F. Lamarque et al. 


Title Page 

Abstract I Introduction 
Conclusions I References 
Tables I Figures 



Printer-friendly Version 
Interactive Discussion 




















ACCMIP, there are considerable differences between the estimated change from 1980 
to 2000 and the simulated one. This is particularly the case over the United States, 
where changes in anthropogenic NH X (and to a lesser extent NO x ) emissions are much 
smaller than over Europe, and do not lead to the observed deposition rate changes. 

In terms of ice-core analysis, the ACCMIP multi-model mean captures many of the re- 
gional features of deposition, but there seems to be an overall estimation in 2000, in 
the Northern Hemisphere Arctic, while Antarctic deposition fluxes are well represented. 

The discussion of the total deposition (wet + dry) in the ACCMIP multi-model mean 
confirms to a large qualitative extent the single-model results discussed in Lamarque 
et al. (2011). In particular, there are large regional increases in 2100 N deposition 
in Latin America, Africa and parts of Asia under some of the scenarios considered. 
Increases in South Asia are especially large, and are seen in all scenarios, with 2100 
values more than double 2000 in some scenarios and reaching region averaged values 
of > 1 300 mg(N) rrf 2 yr _1 in RCP 2.6 and 8.5, ~ 30-50 % larger than the values in any 
region currently (2000). The multi-model mean deposition fields as discussed in this 
study are available as http://acd.ucar.edu/~lamar/ACCMIP/Deposition/all_fields.tar.gz. 

The analysis presented here shows that there is strong potential in using nitrogen 
and sulfate deposition observational data to identify gaps in our understanding of the 
respective precursor emissions. Expansion of the WMO wet deposition assessment 
and more ice cores will provide a window on nitrogen and sulfate emissions, possibly 
helping to understand difficulties for models to reproduce long-term trends such as 
surface ozone in the Northern Hemisphere (Lamarque et al., 2010). 

Supplementary material related to this article is available online at: 
http://www.atmos-chem-phys-discuss.net/1 3/6247/201 3/ 
acpd-1 3-6247-201 3-supplement.pdf. 
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Table 1. List and principal characteristics of ACCMIP simulations (P indicates the primary sim- 
ulations, O the optional ones). SSTs stands for sea-surface temperatures and GHGs for green- 
house gases. Adapted from Lamarque et al. (2013). 
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Table 2. Summary of global totals (Tg(N)yr -1 ) for deposition and emissions in 2000 related to 
the nitrogen oxide emissions. Total dep. (4th column) is the sum of wet and dry deposition while 
Total emi. (6th column) is the sum of NO x surface and aircraft emissions (eminox) and lightning 
emissions (emilnox). Note that deposition is slightly larger than emissions due to the net input 
of nitrogen (usually in the form of nitric acid) from the stratosphere (approx. 1 Tg(N)yr _1 , except 
for the GISS models which have a 5 Tg(N) yr _1 input from the stratosphere). 
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Model Dry Wet Total dep. eminox emilnox Total emi. 

Tg(N)yr _1 


CESM-CAM-superfast 

17 

29 

46 

42 

4 

46 

CMAM 

27 

23 

50 

47 
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51 

GEOSCCM 

12 

33 

45 

40 
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45 

GISS-E2-R 

14 

39 

53 

41 

8 

49 

GISS-E2-TOMAS 

17 

37 

54 

41 
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49 

MOCAGE 
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27 

47 

43 
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48 

NCAR-CAM3.5 

20 

29 

49 

43 
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47 

STOC-HadAM3 
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27 
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UM-CAM 
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Multi-model mean 
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Table 3. Summary of global totals (Tg(N)yr -1 ) for year 2000 of deposition and emissions related 
to the ammonia emissions. Total (6th column) is the sum of wet and dry deposition for NH 3 and 
NH 4 Eminh3 includes anthropogenic and natural (soils and oceans) emissions. 
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Dry NH 4 Wet NH 4 
Tg(N)yr 1 

Total 

eminh3 

GISS-E2-R 

22 

5 

4 

29 

58 

57 

NCAR-CAM3.5 
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Multi-model mean 
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Table 4. Summary of global totals (Tg(S)yr -1 ) for year 2000 of deposition and emissions related 
to the sulfur emissions. Total dep. (4th column) is the sum of wet and dry deposition (both sums 
of S0 2 and SO^ - ) while Total emi. (8th column) is the sum of S0 2 , SO^" and DMS emissions. 
Note that the deposition total should be smaller than emission total since the formation of sulfate 
from the oxidation of DMS has a yield smaller than 1 ; this was estimated to be 87 % in Boucher 
et al. (2003) when both gas-phase and aqueous-phase reactions are taken into account. 
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Wet 

Total dep. 
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emidms 

Total emi. 

CESM-CAM-superfast 
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64 
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19 

83 

GISS-E2-R 

51 

45 

95 

65 
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28 

95 

NCAR-CAM3.5 

24 

55 

79 
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19 

83 

NCAR-CAM5.1 
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56 
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STOC-HadAM3 

41 

43 

84 

69 

0 

20 
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Multi-model mean 
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Table 5. Summary of statistical analysis of the evaluation of nitrate wet deposition (wetno3), 
ammonium wet deposition (wetnh4) and sulfate wet deposition (wetso4) for PhotoComp, HTAP 
and ACCMIP MMM over the 3 main analysis regions. 







wetno3 






North America 



Europe 



Asia 



PhotoComp 

HTAP 

ACCMIP 

PhotoComp 

HTAP 

ACCMIP 

PhotoComp 

HTAP 

ACCMIP 

Linear fit slope 

1.0 

1.0 

0.9 

0.3 

0.3 

0.3 

0.5 

0.5 

0.4 

Linear fit intercept 

-15.4 

-10.3 

-18.7 

16.0 

25.0 

33.2 

22.2 

23.1 

21.3 

Mean bias 

-0.4 

-0.2 

-0.5 

0.4 

0.6 

0.7 

0.4 

0.5 

0.4 

Mean observations 

1.9 

1.9 

1.9 

3.0 

3.0 

3.0 

2.6 

2.6 

2.6 

Mean model 

2.3 

2.1 

2.4 

2.6 

2.4 

2.3 

2.2 

2.1 

2.2 

Correlation coefficient 

0.8 

0.9 

0.9 

0.6 

0.6 

0.6 

0.8 

0.8 

0.8 

Fraction within ±50% 

77.0 

84.3 

68.3 

75.0 

85.2 

85.2 

84.0 

84.0 

88.0 






wetnh4 






North America 



Europe 



Asia 



PhotoComp 

HTAP 

ACCMIP 

PhotoComp 

HTAP 

ACCMIP 

PhotoComp 

HTAP 

ACCMIP 

Linear fit slope 

0.8 

0.9 

0.4 

0.4 

0.4 

0.3 

0.8 

0.7 

0.1 

Linear fit intercept 

-3.3 

-6.3 

17.1 

7.7 

17.3 

46.3 

21.1 

18.8 

77.9 

Mean bias 

-0.1 

-0.1 

0.2 

0.3 

0.5 

0.9 

0.7 

0.6 

1.7 

Mean observations 

1.6 

1.6 

1.6 

3.4 

3.4 

3.4 

4.0 

4.0 

4.0 

Mean model 

1.7 

1.7 

1.4 

3.1 

2.9 

2.3 

3.3 

3.4 

2.3 

Correlation coefficient 

0.9 

0.9 

0.8 

0.6 

0.6 

0.6 

0.8 

0.8 

0.2 

Fraction within ±50% 

82.2 

84.8 

73.0 

73.9 

79.5 

76.1 

76.0 

68.0 

48.0 






wetso4 






North America 



Europe 



Asia 



PhotoComp 

HTAP 

ACCMIP 

PhotoComp 

HTAP 

ACCMIP 

PhotoComp 

HTAP 

ACCMIP 

Linear fit slope 

0.9 

1.0 

0.7 

0.4 

0.6 

0.3 

0.4 

0.5 

0.4 

Linear fit intercept 

-13.0 

-13.9 

-2.4 

19.9 

-11.3 

32.7 

46.8 

36.1 

61.7 

Mean bias 

-0.5 

-0.5 

-0.1 

0.6 

-0.6 

1.0 

2.2 

1.9 

2.7 

Mean observations 

3.1 

3.1 

3.1 

4.0 

4.0 

4.0 

6.9 

6.9 

6.9 

Mean model 

3.6 

3.6 

3.2 

3.4 

4.6 

3.0 

4.7 

5.0 

4.2 

Correlation coefficient 

0.9 

0.9 

0.9 

0.6 

0.6 

0.7 

0.9 

0.9 

0.9 

Fraction within ±50% 

70.4 

70.0 

68.3 

78.7 

52.8 

82.0 

80.0 

88.0 

72.0 
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Table 6. Relative change 2000 minus 1980 (expressed in %, relative to 1980) for NO y wet 
deposition (wetnoy), NO x emissions (eminox) and total precipitation averaged over the region 
of interest. 


Model 

United States 
wetnoy eminox 

precip 

wetnoy 

Europe 

eminox 
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CMAM 
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-2.0 

-0.2 

-4.9 

-13.2 

4.1 

Multi-model mean 

4.9 

0.9 

0.9 

-8.7 

-11.9 

1.6 


g 

</)’ 

o 

c 

< 7 > 
C/) 

O' 


“0 


Q) 

"O 

CD 


g 

c/)‘ 

o 

c 

C/) 

(/) 

o' 


Tl 

0) 

“O 

CD 


g 

c/)' 

o 

c 

O) 

o' 


T) 

0) 

"O 

CD 


6278 


ACPD 

13, 6247-6294, 2013 


Multi-model mean 
nitrogen and sulfur 
deposition 

J.-F. Lamarque et al. 


Title Page 

Abstract I Introduction 
Conclusions I References 
Tables I Figures 



Printer-friendly Version 
Interactive Discussion 




















Table 7. Average nitrogen (NH x + NO y ) deposition (mg(N)m 2 yr 
defined in Lamarque et al. (2011). 

^ over specific regions as 

Discussion Pap 

Region 

2000 RCP26 2100 

RCP45 2100 RCP85 2100 

Canada 

203 

150 

148 

201 

CD 
— ? 

USA 

619 

416 

412 

550 



Mexico 

423 

464* 

351 

503* 


C. America 
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313* 

242 

287* 

g 

</) 

Brazil 
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397* 

299 

465* 

o 

c 

Rest of S. America 
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369* 

253 

359* 

c n 

C/3 

N. Africa 
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139 
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O 

13 

W. Africa 
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524* 

450* 

573* 

“0 

CD 

E. Africa 

291 

461* 

355* 

570* 

“D 

CD 

S. Africa 

295 

389* 

258 

458* 

— * 

W. Europe 

776 

448 

467 

662 

— 

C. Europe 

1054 

685 

577 

1256* 

n 

Turkey 

697 

541 

512 

657 

c/5* 

Ukraine 

852 

584 

415 

1049* 

c 

Kazakhstan area 

294 

361* 

201 

328* 

C/3 

Russia 

287 
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304* 


Middle East 

267 
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~0 

CD 

South Asia 
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1550* 
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1318* 

■a 

CD 

Korea region 
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894 


East Asia 
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Southeast Asia 
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Indonesia 
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Japan 
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c 

C/3 

Oceania 

130 

133* 
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C/3 

o’ 

Greenland 
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39 

46 

58* 

Z3 

“0 

* indicate regions where the deposition increase with respect to 2000. 
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Table 8. Average sulfur deposition (mg(S) m 2 yr 1 ) over specific regions as defined in Lamar- 
que et al. (2011). The deposition does not increase in any of these regions with respect to 
2000. 

Discussion Paper 

Region 

2000 RCP26 2100 

RCP45 2100 

RCP85 2100 

Canada 
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38 
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Fig. la. Nitrate wet deposition (mg(N)m 2 yr 1 ) for 2000. Contours are for the multi-model 
mean, filled circles are for the wet deposition network observations. 
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Fig. 1 b. Ammonium wet deposition (mg(N) m 2 yr 1 ) for 2000. Contours are for the multi-model 
mean, filled circles are for the wet deposition network observations. 
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Fig. 1c. Sulfate wet deposition (mg(S) m" 2 yr _1 ) for 2000. Contours are for the multi-model 
mean, filled circles are for the wet deposition network observations. 
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PhotoComp Asia 





Fig. 2a. Scatterplot of nitrate wet deposition (mg(N) m“ 2 yr _1 ) over North America (left column), 
Europe (middle column) and East Asia (right column). Top row shows the model results from 
PhotoComp, middle row from HTAP and bottom row from this study. See text for details. 
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Fig. 2b. Scatterplot of ammonium wet deposition (mg(N) m “ 2 yr -1 ) over North America (left 
column), Europe (middle column) and East Asia (right column). Top row shows the model 
results from PhotoComp, middle row from HTAP and bottom row from this study. See text for 
details. 
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Fig. 2c. Scatterplot of sulfate wet deposition (mg(S) m 2 yr 1 ) over North America (left column), 
Europe (middle column) and East Asia (right column). Top row shows the model results from 
PhotoComp, middle row from HTAP and bottom row from this study. See text for details. 
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Fig. 3. Change in deposition (2000-1980) over the NADP (left column) and EMEP (right 
column) networks. Top row is for nitrate (kg(N) ha -1 yr _1 ), middle row is for ammonium 
(kg(N) ha -1 yr _1 ), bottom row is for sulfate (kg(S) ha -1 yr _1 ). 
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Fig. 4a. Nitrate deposition from ice-cores (mg(N)m" 2 yr _1 , filled circles, see Table S5) and 
MMM for 1850 (left column) and 2000 (right column). Note different color scales. Top panel is 
for the Arctic region, bottom panel is for the Antarctic region. 
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Fig. 4b. Ammonium deposition from ice-cores (mg(N) m“ 2 yr _1 , filled circles, see Table S5) and 
MMM for 1850 (left column) and 2000 (right column). Note different color scales. Top panel is 
for the Arctic region, bottom panel is for the Antarctic region. 
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Fig. 4c. Sulfate deposition from ice-cores (mg(S) nrT 2 yr , filled circles, see Table S5) and 
MMM for 1850 (left column) and 2000 (right column). Note different color scales. Top panel is 
for the Arctic region, bottom panel is for the Antarctic region. 
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Fig. 5a. Total (wet + dry) NO y deposition 1850-2100 (mg(N) m“ 2 yr _1 ). Top row shows 1850, 
1980 and 2000. Middle row shows 2050 for RCP2.6, RCP4.5 and RCP8.5. Bottom row shows 
2100 for RCP2.6, RCP4.5 and RCP8.5. 
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Fig. 5b. Total (wet + dry) NH X deposition 1850-2100 (mg(N)m 2 yr 1 ). Top row shows 1850, 
1980 and 2000. Middle row shows 2050 for RCP2.6, RCP4.5 and RCP8.5. Bottom row shows 
2100 for RCP2.6, RCP4.5 and RCP8.5. 
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Fig. 5c. Total (wet + dry) SO x deposition 1850-2100 (mg(S) m“ 2 yr _1 ). Top row shows 1850, 
1980 and 2000. Middle row shows 2050 for RCP2.6, RCP4.5 and RCP8.5. Bottom row shows 
2100 for RCP2.6, RCP4.5 and RCP8.5. 
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Fig. 6. Time evolution of deposition acidity (10 6 molm 2 yr ^ computed as 2 SO x + NO y -NH x . 
Blue regions indicate areas where deposition is more basic than it would have been without 
those respective emissions. 
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